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We discuss tiie baryonic contribution to QCD thermodynamics near the QCD phase transition, 
which we split into "stringy excitations" and the "chiral" lowest states. Our finite-T string model for 
the former component is inspired by the lattice data on static potentials and ideas of string survival 
even above Tc. it is used to explain two sets of baryonic susceptibilities calculated on the lattice by 
the Bielefeld-BNL group. Two new ingredients of the model are (i) the near-Tc tightening of the 
strings, and (ii) physical upper limit of the string length. Then we proceed to more subtle effects 
related to chiral restoration dynamics near T^. we suggest that "melting" of the sigma terms of 
the lowest nucleon/A masses can explain these susceptibilities. In a discussion, we consider bound 
monopole-quark states, as a possible explanation to some deficits of contribution at T = (1 — 1.4)rc. 

PACS numbers: 



I. INTRODUCTION 

The issues discussed in this paper deal specifically with 
the sector of the QCD states with a nonzero baryonic 
charge, the baryons in the confining phase and quarks 
in the deconfined one. They are however part of a more 
general discussion about the fate of hadronic states near 
and above the deconfinement phase transition. 

Hadronic spectroscopy studies the lowest hadronic 
states in great details: but as the excitation energy grows 
the task of exact quantum number identification gets 
more challenging. Eventually the decay widths of the 
hadrons become comparable to distances between the 
states, or exceeds it F > \Ei — Ej\: here identification 
of the individual states becomes impossible. Yet it does 
not mean that the states (and their density) do not exist: 
multiple examples from other fields of physics - e.g. the 
overlapping neutron resonances in nuclear physics - show 
that. The thermodynamical quantities are averaged over 
all states excited at some well defined conditions, and 
those include the states which we cannot identify indi- 
vidually. We need certain models about hadronic states 
at one hand, to understand/reproduce the available (lat- 
tice) results on the thermodynamics: this is the main 
logical direction of this paper. 

While the temperature and the baryonic chemical po- 
tential are small T,fi <^ Tc, hadronic matter is a di- 
lute gas of basic hadrons, e.g. pious and nucleons, and 
thermodynamics must be given by elementary ideal gas 
formulae, plus interaction corrections given by the stan- 
dard virial expansion. As one approaches Tc from below, 
many excited states of basic hadrons come into play, and 
hadrons become more densely packed and more strongly 
interacting. It is widely believed that a consistent treat- 
ment of this matter is given by a "resonance gas" approxi- 
mation, in which interaction is taken care of by counting 
all resonances on equal footing with stable hadrons [H. 
Furthermore, when the vicinity of the critical region is 
reached T ^ Tc — 170 — 190 MeT^, new phenomenon (re- 
ferred to as "rise" in the title) happens, in which ther- 



modynamical quantities grow by about an order of mag- 
nitude in a rather narrow interval of temperatures. 

Statistical bootstrap model (hadrons are bags made 
of other bags) lead Hagedorn in 1960's to the conclu- 
sion that the density of hadronic states p(m) may grow 
exponentially with the mass m. If this be the case, Hage- 
dorn argued, at some temperature this growth would be 
cancelling the Boltzmann exponent, preventing hadronic 
systems from being heated to temperatures above it. 
Thus he predicted certain maximal possible temperature 
in Nature, T < Tnagedom [2|- The development of QCD 
in 1973 has lead to the opposite conclusion: the high- 
T phase of QCD is in fact the so called quark-gluon 
plasma (QGP) 3] in which the density of states grows like 
the power of the mass, like for other "normal" systems. 
Later, Hagedorn reinterpreted his argument as the maxi- 
mal possible temperature of the hadronic phase, taking it 
as an early indication to the existence of the QCD decon- 
finement phase transition. The basic physics of Hagedorn 
phenomena was explained by Polyakov and Susskind in 
1970's who pointed out that exponential growth of the 
density of states is the natural consequence of the QCD 
string model. 

As T gets higher than Tc, hadronic (e.g. baryonic) con- 
tribution to thermodynamics gets replaced by free quarks 
and gluons.. The issues we will be discussing in this paper 
is what exactly happens in the so called "RHIC domain" , 
T = {1 ~ 2)Tc. As we will see, lattice susceptibilities tell 
us that the rapid rise of hadronic states changes to nearly 
equally rapid fall, at T ^ l.lTc: there is no consensus 
about its exact physical nature. This question will in 
turn lead us to the issue of whether the very existence 
of strings above Tc is or is not possible. Our point of 
view is that not only it is, but recent RHIC discovery of 
the so called "soft ridge" shows existence of quite robust 
strings, surviving for several fm/c under such conditions, 
for related discussion see We also would claim that 
"stringy" excitations dominate the rise and the fall of the 
baryonic excitations: the latter related to the increase of 
quark density and shortening of the mean length of the 
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strings. 

In order to discuss these issues and provide some 
(model dependent) answers to those and many similar 
questions we had selected a particular subset of hadrons, 
namely the baryons. The reason for that is that by con- 
sidering susceptibilities - the quantities obtained by cer- 
tain number of derivatives over the quark chemical po- 
tential II - we insulate the issue from much more compli- 
cated sector of states with zero baryon number, including 
mesons, glueballs or even more exotic quasiparticles such 
as monopoles/dyons. 

The inputs to our study is provided by the lattice data 
sets from two papers, by UK-Bielefeld-BNL group, to be 
referred below as the Set 1 and Set 2, respectively. The 
former one [5| is for Nf — 2, two quark flavors with a 
rather heavy pion mass, « 770 MeV. The latter 
is a more recent work [6| is closer to the real world, it 
includes the strange quark and also has much smaller 
pion mass « 220 MeV. (Unfortunately, the latter 
one has also somewhat larger fluctuations and thus error 
bars.) Since those two sets correspond to two different 
worlds, all quantities in question ~ including the baryon 
masses and the critical temperatures themselves - are 
different, and reproduction of their differences would be 
one of the tasks for our model. 

The most prominent feature of these susceptibilities 
are illustrated in Fig HI which shows the ratio of the 4-th 
to the second derivative over quark chemical potential, 
to be called kurtosis for brevity. One can see from this 
figure that the kurtosis drops rapidly from 9 to 1 as the 
temperature crosses Tc. It finds a very simple expla- 
nation: while the baryonic partition function is propor- 
tional to cosh(3/x/T), that of quarks has just cosh(/i/T), 
and two derivatives of those functions explain 9 and I in 
the corresponding phases. This means that the confined 
phase is dominated by baryons and the deconfined one 
by quarks. 

While this conclusion is basically correct, a detailed 
study of these transition and the T-dependence of the 
susceptibilities will reveal many more interesting details 
of the transition. The susceptibility in particular show 
a very dramatic rise at Tc, complemented by a rapid drop 
right after that. Both these features are driven by the 
baryon sector, and not by the quarks - as we emphasized 
by the title of the paper. The reason for this to happen, 
as well as some model describing it, is the main part of 
this paper. 

Physics related to these data has been discussed by 
Liao and Shuryak fl\ , who emphasized the presence of at 
least the lowest baryons in the region of T = (I — lA)Tc, 
and in particular the T, /i-dependence of their masses. 
It was based on a previous paper by the same au- 
thors 8] in which lattice-based potentials between quarks 
were used to study when quark quasiparticles can form 
bound states. From a variational solution to the 3-quark 
Schrodinger equation they concluded that only the lowest 
states should survive, remaining bound in the T interval 
just mentioned. Another issue is whether the baryon and 



quark effective masses are /i-dependent. If they are, that 
generates deviations from 9 and I in the kurtosis, and 
certain changes in the susceptibilities: as we will show 
below, those compare well with the latest lattice data. In 
our discussion below we will separate those lowest states 
from "stringy" excitations. 

Weise and collaborators [9l-[ll| have provided very nice 
descriptions of these data in the PNJL model: their 
main idea was that the single quark and diquark con- 
tributions are both suppressed close to Tc by a small 
VEV of the Polyakov loop {L{T)). Ahhough this model 
has no baryons per se, only the colorless combinations of 
three quarks (so to say, "proto-baryons" ) are the color- 
less states which are not suppressed by the Polyakov loop 
in the near- Tc region, and are thus the main contributor 
to thermodynamics there. 

Finally, there are two more important physics issues, 
both related to the chiral dynamics. The first is the 
dependence of baryonic susceptibilities (and of course 
masses) on the value of the light quark mass niq in the 
theory. It has been studied a lot by the lattice com- 
munity, but in the vacuum (at T = 0) rather than at 
T sa Tc. The second is the derivative of the baryon mass 
over the quark mass - known as the "sigma term" . It is 
a chirally odd quantity, expected to vanish at T > Tc in 
the chiral restoration transition. So, we will be look for 
observable consequences of that, in the lattice data set 
under consideration. 

The paper is structured as follows. Section |ll] is de- 
voted to "stringy excitations" of baryons, in it we develop 
a somewhat modernize model of the phenomenon, which 
goes back to Polyakov-Susskind excitation of strings. We 
will show that this phenomenon can indeed naturally de- 
scribe very rapid rise in baryonic contribution to ther- 
modynamics near Tc. We have also added a restriction 
on the upper limit of the string length related to quark 
density: this turns out to produce the observed "fall" of 
the stringy baryons. The density of quarks in the near- 
Tc region is calculated with recent lattice data on quark 
effective masses as well as Polyakov loop, see section ITTll 
In section IIVI we subtract the quark contribution from 
total susceptibilities and compare it to our model calcu- 
lation for "stringy" and "non-stringy" baryons. Effects 
of chiral dynamics, which determine ruq, T and /i de- 
pendencies of the baryon mass, are discussed in section 
IVl Discussion of some open issues is placed in section 
IVII one of them is the role of bound states between 
quarks and magnetic objects, monopoles/dyons, which 
(to our knowledge) are discussed in the QGP context 
for the first time. We present some crude estimates of 
their possible contribution, and find that it seems to be 
of the same magnitude as the "unaccountable" contribu- 
tion with quarks and baryons. 
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FIG. 1: (Color online)The temperature dependence of the 
"kurtosis", the ratio of susceptibihties d^jd^. The Set 1 (2 
flavor QCD at = 770MeV) are shown by black circles, 
while Set 2 (2+1 flavor QCD at = 220MeV) is shown by 
(brown) solid boxes for Nt = 6 and open boxes for Nt — 4). 
The upper and lower horizontal lines correspond to the values 
9 and 1, respectively. 



II. STRINGY BARYONS 

A. Hagedorn transition and the string model 

Specific realization of Hagedorn's ideas came with the 
advent of the QCD strings. Let us briefly remind this 
weh known argument here, for introduction of notations. 
An important parameter of the string a, to be referred 
to as "the string scale", is loosely speaking a length of 
the string section which can be moved independent from 
other segments. (This loose idea implies that it should 
be some 0(1) number times the string width.) Its precise 
definition is given in terms of the number of string states. 
Suppose a string has certain length L and tension a, so 
its energy is E ^ aL. 



which diverges as the temperature approaches the Hage- 
dorn critical temperature 



Tc — aa/ln{5) 



(3) 



Note that at this temperature the exponent (the free en- 
ergy of the string F = E — TS) goes to zero, allowing the 
string to get arbitrarily long. 

Now, does this simple string model describe properties 
of the deconfinement transition observed on the lattice, 
for various gauge theories? 

Before answering this question, one needs to specify 
the string scale scaling. The string tension scales as 



E^a" 



(4) 



if E is the electric field strength and a is the string width. 
At the same time the electric flux, which should match 
the electric charge at the end of the string, is flux — 
Ea? = const, from which one concludes that the scaling 
must be 



E^ l/a^ 



(5) 



Using this argument one expects that Tc/^/a should be 
some universal constants: and in fact for pure gauge the- 
ories with different number of colors Nc this prediction 
actually works well for large Nc, see e.g. [l2| . 

For readers who are not simply satisfied by the scaling 
of Tc and are unsure if the deconfinement is the string 
Hagedorn transition, it is instructive to look at lattice 
measurements of the string free energy. Such measure- 
ments, for Nc as large as 12 has been made by Bringoltz 
and Teper [l^ : see one of their figures is reproduced in 
Fig 12] As one can see from it, with increasing T the ef- 
fective (free energy) tension is indeed decreasing, roughly 
linearly: this would agree with the idea that the energy of 
the string is hardly affected by the temperature and the 
decrease is due to T times the entropy of string states. 
However at T /Tc = 1 the free energy tension is not zero 
but only about half of the vacuum value: using metasta- 
bility of the 1-st order transition the authors were able to 
extend the curve about 5% above Tc and show that the 
trend continues: the projected Hagedorn temperature is 
perhaps Th = l.llTc or so. Thus, the deconfinement 
at large Nc is not the Hagedorn transition, although it 
scales similarly and is numerically close to it. 



p{L) ~ exp{S{L)) ^ 5^/-^ (1) 

Here 5=2d-l is the number of directions the string may 
turn to in d-dimensional space without self-crossing, after 
a segment of length equal to a. Now, combining the 
entropy and energy together one gets stringy partition 
function in the generic form 



B. Strings at finite temperatures in QCD-like 
theories with quarks 

QCD with the physical quark masses is different from 
large Nc pure gauge theories because instead of the first 
order transition with a large jump it has instead rather 
continuous "crossover" transition. The continuity of the 
transition allows us to ask if any role is played in the 
QCD phase transition by the Polyakov-Susskind string 
argument outlined above, at and above Tc- 
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FIG. 2: (Color online) Temperature dependence of the free 
energy tension op, for the SU(IO) gauge theory according to 
Ref.[i2| 



for tensions 

cry = <tf + Tas (10) 
as = -d(7F/dT (11) 

The free energy string tension ap, obtained by solving 
(ITU|) and pT|) with the string tension of the potential 
energy as the input: 

a.^-rjlj-^ (12) 

At To, the free energy tension vanishes, so it is the 
Hagedorn temperature of QCD, which will be fixed later. 
In Fig. ([3]) we show the input string tension of the poten- 
tial energy and the resultant free energy string tension. 
The input are lattice cry extracted from lattice work (T6| 
parametrized by: 

^=26.21e(-"-^°V(^/^c-i)^+6.78xio-) (^3^ 



Before we start answering this question, let us first 
review what is known about QCD strings, at zero and 
finite-T, microscopically . At T = there are extensive 
lattice studies of string tension, shape, internal structure 
[Tsj and even the string breaking It is done by in- 
serting two static charges and studying their influence on 
the gauge fields: it helps by the fact that at T = the 
string is just straight, with small influence of quantum 
vibrations. We will not go into details here, just men- 
tion that the shape fits well into the "dual superconduc- 
tor" picture, with condensed monopoles creating a coil 
of magnetic current around it. The distribution of the 
electric field fits well into the solution of (dual) Maxwell 
equations in superconductor found by Abrikosov 

E{r) ^ Koir/X) (6) 

with the length parameter A w .15 fm, {1.3GeV)~^ 
[isj . Unfortunately, at finite T thermal fluctuations of 
the string make such studies much more challenging, and 
so we do not have similar results in the near-Tc region. 

Yet there are still extensive lattice studies of the ther- 
modynamical quantities associated with a pair of two sta- 
tionary qq [Tsl Il6j. We will not review this large subject 
here, but only comment on the linear (in distance) part 
of these static potentials. Let us define the free energy 
tension ap{T) and entropy tension ug{T) as follows 



V^av{T)L + Vo (7) 
F ^aF{T)L + Vo (8) 
S = as{T)L (9) 
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FIG. 3; (Color online) The temperature dependence of two 
string tensions, ^^(red dashed line), (TF(blue solid line). The 
entropy term Tas is shown by the (green dash-dotted line). 
The lattice data points shown are for av extracted from lat- 
tice work ^16|]. The constant used in the plot is To — 0.95Tc. 
Its determination will be explained later. 



C. T-dependent string model and density of states 



where Vq is a constant. Standard thermodynamical rela- 
tions U = F+TS, S = —dF/dT imply the same relations 



One can build a simple model based on the same num- 
ber of state as for the fundamental string. Suppose T- 
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dependent string scale a(T) be known: then we postulate 
the density of states to be 



p{L) ^ S^/'^f^) = e^-(^)^ 



(14) 



but now with effective T-dependent string scale. To 
write the density of states per unit mass, we note m = 
mo + av{T)L{T), with mo the contribution from the 
quarks at he ends of the string. As the string tension 
cry(T) changes with temperature, the string length L{T) 
corresponds to certain baryon state changes qualitatively 
in the opposite way, i.e. a tight string is short in length. 
As a result, the quark contribution will also change with 
temperature. In order to take into account this effect, we 
introduce a factor x in the exponent. We further assume 
it is a constant in the temperature interval we explore: 



p(m)|T=o - ^^'^^^ 
p("^)|T>o = ce 



oO"s/(Tv/(m-mo) 



ce ■= = ce 

xas(T)/av(T){ra-mo) 



(15) 



The normalization constant c in ()15|) will be fixed by 
matching the experimental data at T = 0. 

In order to make comparisons with the output of lat- 
tice simulation, where large unphysical quark masses are 
used. We adopt the following parametrization for the 
baryon masses: 



ms = av{T)L{T) + iV,m, + kml 



(16) 



in which the first string term is complemented by 
two more terms, which describe the dependence of the 
"string's ends" on the strange and light quark masses. 

The second term is just additive dependence on the 
strange quark mass, with Ns the number of strange 
quarks in the baryon. The light quark mass (pion mass) 
dependence is included via the slope 



k = 



dniB 



(17) 



known also as a "sigma term". It depends nontriv- 
iall y o n the baryon in question, see lattice studies such 
as [13] ■ The N and A channels have stronger depen- 
dence than the A and S channels, which contain more 
strange quarks. We list the slope k extracted from [l3| 
for relevant channel in TablH The fl channel contains 
only strange quarks, therefore we assume the slope to 
be zero. Note the lowest states of each channel is sepa- 
rated from the other resonances by a significant gap, and 
they have different properties (we will call them "chiral 
baryons" and further justify this separation in sectionlV]). 
These lowest states (for three flavors the 56-plet, spin 1/2 
octet Af(938),A(1116), S(1195),S(1317) and 3/2 decuplet 
A(1232),S*(1385),S*(1530),r2(1672), and for two flavors, 
the 20-plet of TV and A) are treated separately from the 
rest of "stringy" (also called hybrid) excitations. 



TABLE I: slope k for different channels 
Af(GeV-i) A(GeV"-^) A{GeV-'^) E(GeV^-i) E{GeV- 



0.79 



0.73 



0.61 



0.55 



0.3 



100- 
90- 
80- 
70- 

p(m)r. 60- 




Pi'n)T 



100- 
90- 
80- 
70- 




FIG. 4: (Color online) The "string inspired" exponential fit to 
the baryonic density of states, for Set 1 and Set 2 (upper and 
lower figure, respectively). The (red) dashed line is the sum 
of over all resonances listed in Particle Data Book ^T§\ , with 
masses corrected to lattice quark mass by expression (|16|l . The 
(blue) solid line is the exponential fit: p(m) = O.lle'"''^'^ (Set 
1) and p(m) = 0.35e'"/^ '^ (Set 2), with unit set properly by 



Let us return to the density of states of "stringy exci- 
tations" , by which we count all resonances to which, for 
smoothness we assign a width of 200A'IeV. We modify 
the baryonic resonances according to (jl6p and fit the den- 
sity of states by the string inspired exponential function 
. Fig|3] shows the exponential fitting to the densities 



6 



of states. 

Comparing our stringy density of states with lattice 
modified baryonic resonances, we have the fohowing 
parametrizations ioi Nj = 2 c = 0.11e™°/^-^(Set 1) and 
Nf = 2 + 1 c = 0.35e'"»/i-^(Set 2). From (HH), we should 
have Too — km^ + Nsiris, which varies in different chan- 
nels. Since we have counted all stringy baryons indis- 
criminately in the density of states for statistical reasons. 
We simply use the smallest toq among all channels. Fur- 
thermore the Hagedorn temperature Tq in will be fixed 
from the relation crs/ay = 1/1. 5T^^^ at T = 0, giving 
To = 0.95Tc. It is important to note the critical temper- 
ature corresponding to Set 1(2 flavors at TOtt = 770MeV) 
and Set 2{2+l flavors at to^ = 220MeV) are slightly dif- 
ferent. The former has Tc = 220MeV, while the latter 
has Tc = 200Afey(They can be obtained fromfHi^). 
Special attention should be paid in converting the abso- 
lute unit to units of Tr, different for each set. 



III. QUARKS IN THE NEAR-Tc REGION 

We postulate that the free energy (pressure) due to 
quarks and antiquaries is given by the usual free gas for- 
mula 

p,(T,/i) = d^mlTK2{m,/T)cosh{p/T){L{T)) (18) 

modified by the last factor, the so called the Polyakov 
loop. 



1 



-Tr 



(19) 



induced by the presence of nonzero Aq gauge fields. The 
corresponding lattice data and our parametrization of 
both the assumed quark effective mass and the Polyakov 
factor are described in the Appendix. 

As those calculations are rather standard, we would 
not show the quark contribution itself, showing in several 
plots below only the susceptibilities with and without 
quark contribution, shown by boxes and circles, respec- 
tively. The latter contributions show a near- Tc peak and 
then decrease, disappearing around 1.5Tc. Our task in 
the section to follow would be the interpretation of these 
"non-free-quark" contributions. 

As this and similar examples show, one has inevitably 
apiT) crossing zero. We would however argue that un- 
like in the string-Hagedorn transition, the string length 
should not be allowed to become indefinitely large and 
the partition function to diverge. The reason for that is 
string breaking/flipping. Let us introduce the quark scale 



-.,-1/3 



(20) 



related to their density in the QCD matter at tempera- 
ture T. 

Note the quark density Uq from naive derivative of the 
quark pressure with respect to /i will give sinh(/i/T), 



which vanishes at /i = due to the cancellation between 
quarks and antiquarks. We instead use cosh(/x/T), as 
both quark and antiquark can be the end points of a 
string. It is also required by the —fi symmetry. 

The degeneracy factor d = 2(25 -I- l)NcNf. Suppose we 
start with a string, whose ends are fixed, for simplicity 
at some nearby points. If the string length is L and its 
shape is random, the average distance from the original 
ends would be given by random walk expression 



< {6xf >- a^{L/a) 



(21) 



where L/a is the number of random turns. We will then 
argue that as a string reaches the "domain of influence" 
of other quarks, the string would not grow anymore in 
length but get switched to another quark. We thus pos- 
tulate that 

<{Sxf >^a\L/a)<Ll (22) 
sets the upper limit in the integral over the string length: 



Lmax — L^q/o, 



(23) 



with y is some constant of order 1. Correspondingly, 
mmax = 'mo + o'v{T)Lmax- The temperature dependence 
of Lq is such that at T < Tc it is very large, reaching large 
but finite value at Tc, and then rapidly decreasing. 

We compare various length scale appearing in the 
model in FiglS] At temperature T ~ 0.8Tc, the string 
length can get very long, but long strings are suppressed 
by Boltzmann factor. As the temperature crosses the 
critical value, the density of states starts to win over the 
Boltzmann factor and the string length cutoff becomes 
important, as illustrated by the dropping of the Lmax 
near Tc. At T ~ l.lTc, the stringy baryon picture breaks 
down as Lmnx ^ a. 



IV. THERMODYNAMICS OF THE "STRINGY 
BARYONS" IN THE NEAR-T^ REGION 

Now we collect all the pieces of our string-inspired 
model and calculate several susceptibilities (i„(T), to be 
compared with Set 1 and 2 lattice data. The expressions 
of the pressure and susceptibilities due to baryons are 
shown in the following: 



Td\iiZ{T,^i) 
dV 

1 

^T2 
' ~ 9(/i/T)" '^=° 



J2±!dK,im,/T)cosh{3f,/T) 



TT- T2 

K2{m/T)cosh{3fi/T)p{m)dm 



(24) 
(25) 



The lowest states counted separately by the index i in- 
clude the 20-plet of N, A for Nf = 2 or the members of 
the 56-plet, if the strange quark is present. The second 
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FIG. 5: (Color online) A comparison of the string length cutoff 
imaa;(blue dashed line) and string scale a(red solid line) at 
different temperatures for two ffavors QCD. (green dash- 
dotted line) is also included for comparison. The parameters 
used in the plot is y = 0.08, To = 0.95Tc 



term is the contribution from the continuum assumed to 
be stringy excitations. The upper bound in integration is 
given by mmax = wo + o'v{T)Lmax and the lower bound 
nT-min IS givcn by the mass of the first excited state in N 
channel, with unphysical quark mass properly taken into 
account by ([TBI) . 

Now is the time to show the results and confront it 
with lattice findings on susceptibilities. We show those, 
with and without the free quark contribution, in Figl6] 
for the first susceptibility d2. 

The red dashed line shows the contribution of the 
"stringy" baryons, while the dash-dotted line shows the 
contribution of the lowest "chiral" states. First of all, the 
stringy states, albeit higher in mass, predominates over 
the lowest states near Tc, and their rise explains the data, 
both in magnitude and in the growth rate. Their domi- 
nance is more prominent at larger bare quark mass(Set 1) 
than at smaller bare quark mass (Set 2), the latter having 
smaller masses for the lowest states. Another reason why 
Set 2 has larger contribution of the non-stringy baryons 
is that it has Nf = 2 + 1 rather than 2, causing further 
proliferation of baryons relative to quarks. 

As the temperature rises above Tc, the string model 
predicts the fall of the baryons caused by quark prolif- 
eration and a decrease of the string length. It happens 
when rrimin ^ "rumax which occurs at T ^ 1.05Tc. This 
value is in qualitative agreement with the result obtained 
from the L^ax a condition in the previous section. Un- 
fortunately, this fall is much more rapid than the data, 
and the question arises what is the physical nature of 
states between 1.1-1.4 Tc for set 1 and 1.1-1.8 Tc for Set 
2. 

Those cannot be the "non-stringy" baryons, shown by 
the dash-dotted curves. First of all, their contribution 
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FIG. 6: (color online) Various contributions to the susceptibil- 
ity d2(r) as a function of the temperature. The upper figure 
shows the lattice data of the Set 1(2 flavors at m.^ = llQMeV) 
from ref. with quark contribution subtracted, are shown by 
(black) circles. The unsubtracted data are also included, as 
(brown) boxes for comparison.) The (red) dashed line cor- 
responds to hybrid (or "stringy") baryons. The (blue) dash- 
dotted line corresponds to the sum of the lowest 20-plet states. 
The (black) solid line is the overall baryonic contribution. 
The parameters used are x = 0.62, y = 0.08. The lower plot 
is for Set 2 (2-1-1 flavors at = 22QMeV) d. The lattice 
data on susceptibility^] with quark contribution subtracted 
are shown as black circles (solid ones for A'^^ = 6 and open 
ones for Nr = 4). The unsubtracted data are also included as 
brown boxes (solid ones for At- — 6 and open ones for At — 4) 
for comparison. The parameters used are x — 0.51, y — 0.08. 



have the wrong shape: it is increasing rather than de- 
creasing. Furthermore, we terminated those lines at the 
melting temperature of baryons calculated in , modi- 
fied for the particular quark masses used here. 



8 



V. CHIRAL DYNAMICS AND THE BARYONS 

A. Why do we have two types of hadrons? 

In a narrow sense, chiral dynamics is related to the de- 
pendence on the dynamical quark masses in the theory 
rriq when they are small or vanishing. In a more broad 
sense the issue includes an old question of to what ex- 
tent the baryons are made of their "pion cloud" , and to 
what extent of quarks and confining strings, as in the 
nonrelativistic quark model. Hadronic physics has been 
oscillating between these two conflicting points of view 
for decades, with drastically different expectations for T 
and fi dependencies of hadronic masses following from 
them. 

One extreme point of view, following from the Skyrme 
model and application of the QCD sum rules, suggested 
that hadronic masses are due to sigma/pion bags and 
are impossible without the chiral symmetry breaking and 
nonzero {ipip)- If so, one would expect drastic changes or 
even vanishing of hadronic masses at the chiral restora- 
tion transitions (2l| . 

The opposite point of view is that the sensitivity of 
baryon masses to chiral dynamics is restricted to explicit 
chirally odd coefficient of the quark mass, the sigma term. 
Here we will adhere to this conservative point of view. 

Before turning to the chiral effect on the baryon mass 
and density of states, we want to point out a significant 
difference between (i) the lowest baryons (in each chan- 
nel, such as iV, A) and (ii) the bulk of hybrid baryons. 

The former ones (i) - as experiment and lattice simu- 
lations told us - are strongly interacting with pions and 
chiral condensate. They do not have nearby parity part- 
ners. Significant part of their mass is believed to be re- 
lated to "the pion cloud" , the extreme model on that is 
of course the Skyrme model which look at them as be- 
ing pionic solitons. Instanton liquid model explains their 
masses well, using propagators based on instanton zero 
modes. 

The latter (ii), the bulk of excited baryons, are on 
the contrary very insensitive to chiral physics. The most 
vivid spectroscopic evidence for that is parity doubling 
phenomenon, see extensive set of data on that in [g^l . As 
was repeatedly emphasized in discussions of the matter, 
close pairs of opposite parity also imply weak coupling 
to pions, which was also evident from the magnitude of 
decay matrix elements. 

Why are these excited baryons behave like this? Our 
interpretation of them as "stringy" excitations would 
suggest that perhaps interaction of the color string and 
chiral condensate is small. This was in fact demonstrated 
on the lattice already 20 years ago [1^ : in the vicinity of 
static quark the (V'V') is reduced by few percents only. 

Instanton liquid model provides another interpretation 
of this phenomenon: it is a reflection of the Dirac spec- 
trum in the field of an instanton, in which there is one 
(chiral) zero mode plus a continuum of non-chiral nonzero 
ones. The lowest hadrons are in this model viewed as 



exceptional states, built from collectivized zero modes, 
unlike all the excited states. 

The lowest excitations are described via some poten- 
tial models (especially accurate for heavy quarkonia) as 
radial and orbital excitations of moving quarks. Excita- 
tions in which quarks are basically at the same state as 
in the lowest hadrons but a string is excited are called 
"hybrid states": there is vast lattice and experimental 
literature about them. The string-Hagedorn argument 
mentioned above implies, that such states are dominant 
at large excitations. 



B. Chiral restoration transition and the lowest 
baryons 

In the model proposed we assumed that the strongest 
sensitivity to small quark mass is entirely due to the "low- 
est" states with each quantum numbers. The mass of 
some baryonic state can be expressed by 



rriB 



{B\il:il)\B) + "gluonic part" 



(26) 



where the coefficient of the quark mass is the scalar quark 
density: it is known as the sigma term. 

In the calculation above we had fixed its magnitude 
from lattice works, in the vacuum: now is the time to 
discuss its dependence on the temperature. Following 
the idea that any chirally-odd quantity should disappear 
at chiral restoration transition, and furthermore that the 
sigma terms for any hadron is more or less proportional to 
the VEV or quark condensate, we propose to incorporate 
the chiral effect and fx dependence through the following 
formula: 



2 ,B\^^\B)iT,^i) 



[B\i,iP\B)\T=^=o 



N,ms + (Tv{T)L{TQn) 



In the absence of knowledge of the expectation value of 
ipil) with respect to the baryonic state, we simply assume 
normalized chiral condensate associated with some bary- 
onic state is given by its counterpart in vacuum state: 



{B\i^m){T,y) (^V)(T,Ai) 



(i3|VV'|i3)|T=M=o (V'V')|t=m=o 



(28) 



The teniperature dependence has been measured on 
the lattice [2J]. We fit the normalized chiral condensate 
by the following: 



{^i^){T,^, = {)) 

(■0-0)|t=m=O 



0.55 - 0.42tanh(9.47— - 9.55) (29) 



The /i dependent term enters quadratically by the 
O — /i symmetry. We will use the collective coor- 
dinate ^(r/rc)2 + (Ai^c)2 instead of T/T^. Similar 
parametrization is also used in 0. We leave 
free parameter. With chiral effect incorporating both T 
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FIG. 7: (color onlme)Tliree figures display the temperature 
dependence of the second, fourth and sixth susceptibilities. 
The points are the Set 1 lattice data with the quark contribu- 
tion subtracted. The lines correspond to the model discussed 
in the text, with /i dependent stringy and chiral baryons. The 
parameters used are x = 0.6, y — 0.08, and /ic = oo(black 
solid line), fic — 2rc(red dashed line) and /ic — 1.33rc(blue 
dash-dotted line). 



and fj, dependence, we are ready to calculate the suscepti- 
bility using (1^ . Since fi appears in many places in 
taking derivatives becomes cumbersome. We simply use 
numerical differentiation method. Further simplification 
can be made by noting the evenness of ([24]) as a function 
of fJ-/T. 2n-th derivative with respect to fi/T amounts to 
n-th derivative with respect to (/i/T)^. 

As already discussed in 0, , in the model of baryon 
gas the deviations of ^4/^2 from 9 follows from the fact 



FIG. 8: (color online) The temperature dependence of the 
second and fourth susceptibilities, for Set 2 lattice data, with 
the quark contribution subtracted. The (black) solid circles 
correspond to A'^t = 6 and open circles to Nt = 4. The curves 
are for the model described in the text, with parameters values 
equal to x = 0.51, y — 0.08, and /ic = oo(black solid line), 
fj,c = 2rc(red dashed line) and /ic — 1.33rc(blue dash-dotted 
line). 



that their masses are effectively dependent on /i. In our 
string model, this appears as a /i dependent density of 
states p{T, fi)^. 

We show the effect of /i dependence of baryon mass 
in FiglTtSct 1) and Fig|Si;Set 2). We see the inclu- 
sion of this effect hardly change the c?2 plot, but shifts 
the peak upwards in and dg plots. Smaller values 
of result in higher peaks. However the peaks of d4 
and dg are over predicted either with(/.tc = 1.33, 2Tc) or 
without(/ic = 00) ^ dependence of the baryon mass. The 
shift of the peak is more prominent near Tc, where the 
chiral condensate drops rapidly, thus the masses of the 
lowest states change rapidly as well. It is also worth not- 
ing that the shift is less prominent in case of 2-|-l flavors 
model. This is because in that case a smaller fraction of 
the mass has changed due to a smaller sigma term. 

It is interesting to note that the peaks of susceptibility 
data(with quark contribution subtracted) move to lower 
temperature as we go to higher cumulants, e.g. ^2 in 
FigEl peaks at T - LIT^ while dg peaks at T - l.OTc. 
This might be an artifact of our simplified treatment of 
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quark contribution (jl8|) . It is worth mentioning that a 
recent study with holographic model f26] shows nontriv- 
ial temperature dependence of susceptibilities of different 
cumulant, which deserves further understanding. 

It is also instructive to plot the kurtosis from our 
stringy baryon model. We show the plot of kurtosis both 
with and without the chiral effect in Figl^l A sharp drop- 
ping of the kurtosis is observed near the critical temper- 
ature. However the limitation of the model prevents us 
to go beyond 1.05Tc, where the kurtosis is supposed to 
drop to 1. The importance of /i dependence of the baryon 
mass is clearly illustrated in FiglHl Without the fi depen- 
dent density of states, the drop is a simple vertical fall, 
while the inclusion of chiral effect produces a bump where 
the chiral condensate drops rapidly. The bump is more 
prominent for 2 flavors Set 1, in which case the lowest 
states have relative large fraction of mass that is chirally 
sensitive. 
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FIG. 9: (color online)The temperature dependence of the kur- 
tosis ^ All discussed lattice data are shown here: Set 1 (2 
flavor at = 770MeV) by (blue) circles, Set 2 (2+1 flavor 
at m-n = 22QMeV) by (red) solid boxes and open boxes, for 
Nt = 6 and Nt = 4 respectively. Four curves in the plot are 
for 2 flavor at = 770MeV without (blue dotted) and with 
(blue dash-dotted) chiral effect, 2-1-1 flavor at m-n = 220MeV' 
without (red solid) and with (red dashed) chiral effect. The 
same parameters are used as in Fig |6l7l8l In particular, we 
have chosen /ic = ITc for cases with chiral effect) 

The sign and magnitude of the chiral restoration ef- 
fect in our model is in sharp contrast with the PQM 
model [l^l, predicting sharp and very large deviation of 
the kurtosis from 9, up to about 15. That happened be- 
cause in Brown-Rho-type dropping of hadronic masses 
prior to chiral restoration transition. In Figl^] there is 
certain mixed signal here: Set 2 data at N-r = 4 and 
Nt — 6 show somewhat different trend. Unfortunately 
current data are still insufficient to make clear distinc- 
tion between our conservative approach (only the sigma 



terms melt at chiral restoration) and the Brown-Rho be- 
havior. We hope this interesting issues will be clarified 
in the future. 



VI. DISCUSSION 
A. Quarks bound to monopoles? 

As we have already emphasized above, lattice data on 
the deconfinement transition give rather different width 
of the transition region for baryons (discussed in the pre- 
ceding section) and the quarks. While the baryons show 
rapid rise and then rapid fall already at about l.OSTc, 
quark suppression is dominated by the Polyakov loop, 
which only goes to 1 at 2Tc. Such mismatch of the widths, 
be nearly and order of magnitude, leave some "hole" in 
the middle: indeed, one can see from Figl7]that our model 
has certain deficit at T = (1.1 — \A)Tc, as compared to 
the lattice data. So, it can be a deficiency of the model, 
or some missing component. 

Thinking about possible missing states leads us to pro- 
pose one idea, which is by itself well known in theoretical 
literature, yet (for our knowledge) never been used in 
QCD thermodynamics context so far. Static monopole 
solution has 3-dimensional fermionic zero mode. It is 
important to note, that its existence follows from topo- 
logical theorems and thus does not depend on details of 
the monopole profile. In particular, it exists whether 
monopoles can or cannot be described semiclassically. 

In vacuum T = this simply means that a quark can 
be bound to the monopole. For a massless quark in vac- 
uum there are two degenerate states, with fermionic num- 
bers 1 and ofs^l- For Nf quark flavors each can be either 
present or absent, leading to 2^^ states. 

Euclidean finite T formalism breaks the degeneracy of 
this situation: the bosonic objects should have periodic 
wave function, and thus time-independence is unaffected, 
while fermionic objects should have antiperiodic wave 
functions in Euclidean time. Thus one should inevitably 
introduce exp{iET) with the lowest odd Matsubara fre- 
quencies E = ±ttT into the fermionic wave function, 
making the composite fermionic magnetic object heav- 
ier 

Mf„, « M™o„o + ttT (30) 

The mass of the monopole is believed to be minimal at 
Tc, while its value is not yet measured. It has been ar- 
gued from BEC condition [2^ that it may be as small as 
Mmono ~ 200 MeV: if so, the minimal mass of the com- 
posite object is about Mfm{T = Tc) ^ .8GeV, which 
is accidentally close to lattice data on quark quasiparti- 
cles. Note however, that such composite has only one 
spin state. 

At temperatures T ^ Tc, the monopoles remain 
strongly correlated and behaves as liquid 29], preventing 
us using ideal gas formula. The scaled monopole density 
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is measured in SU{2) gauge theory on the latticeflOl and 
is parametrized by: 



T3 - /,2(T) l^^J 

with A = 0.557 ^ = 2.69 (32) 
A 

The density of monopole(anti-monopole)-quark bound 
state is estimated as ne~'^ . The additional BoUzmann 
factor arises from the Matsubara frequency. We plot 
the bound state contribution in FigfTUl Taking into ac- 
count a very schematic model and preliminary data on 
the monopole density, we conclude that magnetic com- 
posite fermions do have a potential to be the "missing in- 
gredient" mentioned at the beginning of this section. Of 
course, to know it for sure much more work is needed. In 
particularly, existence of such composite fermions with a 
magnetic charge can be investigated on the lattice, by ob- 
serving the correlations between the monopole paths and 
the fermionic operators (especially the fermionic spin). 
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FIG. 10: The contribution from monopole(antimonopole)- 
zero mode bound state. We have used a degeneracy factor 
2NcNf since the zero mode only has one spin 



B. Outlook 

We hope that multiple details of the last sections have 
not detracted the reader from the main ideas of the 
model, which are (i) exponential rise of the stringy par- 
tition function near Tc, as well as (ii) the upper cutoff of 
the strings length defined by the quark density. Those 



two phenomena are the basis for "the rise and fall of the 
baryons" put into the title. 

Our conservative treatment of the chiral dynamics 
seems to be working. The dependence of the baryon 
masses on m^, light quark bare mass, and related effects 
due to chiral symmetry restoration and vanishing of the 
"pion cloud" are small but interesting corrections to the 
main picture. We have tried ~ using lattice data when- 
ever those are available - to parametrize those effects and 
include them into the model, in the last few sections. The 
reader may at this point ask what is actually achieved by 
doing this. 

Any model is as good as ideas put into it, and in order 
to test it one can see if it has any predictive power. By 
this we mean its application to wide range of the QCD- 
like theories. So far we considered two directions: (i) 
variable light quark mass, and (ii) increasing fi/T or - 
which is basically the same - including the susceptibilities 
with more derivatives. 

We suggest to address the same issues using other theo- 
ries, especially with dynamical adjoint fermions. In this 
case the critical points of the deconfinement transition 
and that of chiral restoration are very different 



decon f inement hiral (33) 

so that there is a new intermediate phase, that of "de- 
confined constituent quarks" . For three colors the decon- 
finement is in this case the first order transition, while 
the chiral one is of the second order, see e.g. [3l[ and 
references for earlier papers therein. Another reason one 
may be interested in this direction is possible relation to 
supersymmetric theories. 

For simplicity and continuity of the observables, one 
may introduce non-dynamical fundamental quarks which 
would have their own chemical potential and form 
baryons with the usual quantum numbers, basically re- 
peating the same analysis as was done for Set.l and 2 we 
used , but with a different ensemble of gauge fields, now 
with some dynamical adjoint fermions. In this case one 
expect to find the "rise and fall" of the stringy baryons 
be seen at the deconfinement transition, while the ef- 
fects related with chiral condensate (and sigma term) 
melting occurring at completely different chiral restora- 
tion temperature. Large intermediate phase would be 
a great place to finally find out whether confinement 
is or is not required for the existence of the "chiral" 
baryons. Another difference between adjoint and funda- 
mental fermions is a different number of fermionic zero 
modes of monopoles/dyons: there are two bound states of 
adjoint fermions rather than one. So, if such states play 
a role above the deconfinement, they should be signifi- 
cantly more important in the adjoint theories, especially 
those with more than one fermion fiavors. 
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Appendix 



In this appendix we collected some information we used 
from different lattice works, to parametrize several quan- 
tities of interest. 

The Polyakov loop in two flavors QCD is shown in 
FigCH It reaches 1 at around 2.5Tc, where the quark 
density reduces to the free quark situation. The Polyakov 
loop is parametrized by: 



{L{T)) = 



2.24 



0.038 + 5.20(r- 0.77) 
-0.164+ 1.048(T- l.Of 



T < 1.07T, 
T > 1.07T, 



=(34) 



The thermal quark mass as a function of the tem- 
perature is shown in Fig ll2l The thermal quark mass 
niq around the critical temperature is important as it af- 
fects the cutoff in the baryon density of states via n,. 
Recent lattice measurement [s^ indicated ruq/T ~ 0.8 in 
temperature range T ~ 1.25 — STc- Thermal mass below 
1.25Tc is still not known. In the absence of additional 
data, we estimate the quark thermal mass as half of the 
internal energy for quark-antiquark at infinity separation 
Uao ■ The data of Uoo are taken from Table I of [l^ [s^ . 
The thermal mass of the quarks are fitted with: 




FIG. 12: (Color online) The points are lattice data on the 
temperature dependence of the effective quark mass, deduced 
from the potentials studied in Ref . [l6|] . The (red) line is a 
fit with expression psp . Direct measurements of the quark 
effective mass from the dispersion relations of the quark quasi- 
particles for two temperatures, T = 1.25rc, l.STc, STc, from 
Ref. [sjl are shown by (blue) circles for comparison. 



m„ 



4.76tanh(12.38|;^ 



13.16) (35) 
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FIG. 11: (Color online) The points are lattice data on the 
temperature dependence of the Polyakov loop from Ref. [isl] . 
fitted with the expression p4l) 



FIG. 13: (Color online)The temperature dependence of the 
chiral condensate, normalized to its vacuum value. The 
(black) dots are taken from lattice work Ref. 24] and fitted 
with the red curve according to expression (1291) . 



The chiral condensate from the lattice measure- 
ment (p9|) . is shown In Fig US] It is defined as, 

-7— 7T ml , 7 , . — and measure at unite mass ratio of 

light and strange quark We will however use this 
value as our normalized chiral condensate, which is valid 
at ^ ^ 

rris 

We fix the dependence by the fohowing ansatz: 



(ViV)(r,M) 



0.55 



0.42 tanh(9.47A/(T/Tc)2 + {ji/iicY - 9.55) (36) 
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